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Abstract 

The recent development of perturbative QCD (PQCD) fragmentation functions 
has strong impact on quarkonium production. I shall summarize B c meson production 
based on these PQCD fragmentation functions, as well as, the highlights of some recent 
activities on applying these PQCD fragmentation functions to explain anomalous J/ip 
and rp' production at the Tevatron. Finally, I discuss a fragmentation model based 
on the PQCD fragmentation functions for heavy quarks fragmenting into heavy-light 
mesons. 



I. Introduction 

One of the biggest ever discrepancies between theoretical predictions and experimental data 
was J ftp and ip' production observed by CDF M at the Tevatron. With the recent devel- 
opment of perturbative QCD (PQCD) fragmentation functions the experimental data can 
be accommodated within reasonable uncertainties. In the following I shall briefly describe 
the essence of PQCD fragmentation functions [@,|],f|], then summarize B c meson production 
based on these PQCD fragmentation functions |||6],[7j and the highlights of some recent ac- 
tivities on applying the PQCD fragmentation functions to explain anomalous J/ip and ip' 



production P,|9| J10| , pdlJl^ , P^ , P^ , p!5| . Finally, I shall discuss a fragmentation model |lq , |17| ,|r 
based on these PQCD fragmentation functions for the fragmentation of heavy quarks into 
heavy-light mesons, e.g., c — > D, D* and b — > B,B*. This model is more attractive than 
previous fragmentation models since it is based on PQCD and the PQCD fragmentation 
functions have the correct heavy quark behavior. 

In general, fragmentation of quarks and gluons lies in the nonperturbative regime so that 
the fragmentation functions cannot be calculated from first principle. But there is a partic- 
ular class of fragmentation functions, namely those for heavy quarks or gluons fragmenting 
into heavy-heavy bound-states, that is calculable in PQCD. Heavy-heavy bound-states refer 
to heavy quarkonia, (cc), (bb), and (be) mesons. To visualize let us consider the hadroniza- 
tion of a heavy quark Q into a meson Qq, which is schematically shown in Fig. [I]. It is the 
lowest order 1-gluon exchange diagram that describes the hadronization process. As shown 
in Fig. [I], the central part is the creation of the quark-pair qq out of the vacuum, followed 
by the binding of q to Q to form the meson Qq. Therefore, the natural scale of this process 
is of order m q , specifically we choose it to be 2m q . Figure |I] could be used to picture the 
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Figure 1: Schematical diagram for a heavy quark Q fragmenting into (Qq) and q. 



fragmentation of a heavy quark into a heavy-light meson when Q = b,c and q = u,d. The 
natural scale is then equal to 2m u> d, which is close to Aqcd, so that we do not expect the 
process to be calculable by PQCD. The nonperturbative physics involved demands a model 
calculation. 

On the other hand, figure [l] can also be used to describe the production of (bb), (cc), 
and (be) mesons when Q = b,c and q = b,c. But the main difference is that the natural 
scale of the process is now of order 2m^ c , which is much larger than Aqcd- in other words, 
the fragmentation process can be reliably calculated as an expansion in a s (2m6 )C ). While 
the diagram in Fig. [I] represents the lowest order a 2 term; higher order corrections can be 
systematically calculated. The arguments for the calculability of gluon fragmentation into 
charmonium and bottomium within PQCD are essentially the same. 

Even though we can calculate the fragmentation process by PQCD, there are bound-state 
effects that have to be taken care of. The nonperturbative bound-state effects can be param- 
eterized by, e.g., wavefunctions or derivatives of wavefunctions of the bound-states at the 
origin. For the case of S-wave mesons there is only one nonperturbative parameter - R(0), 
the wavefunction at the origin; while for P-wave mesons there are two nonperturbative pa- 
rameters, which correspond to the color-singlet and color-octet mechanisms. Fragmentation 
into different spin-orbital states can be obtained by using the appropriate spin projections. 

The b, c quark fragmentation functions can be calculated by this expression 

D ^ Qi(z) = i f dse L-Ml.J±) , im m (1) 

V W\ lg7r 2 J y z 1 - Z J PQ /mQ-*oo \Mq\ 2 

where Q,q = b,c, M = uiq + m q is the mass of the meson, M. is the amplitude for producing 
a Qq and q from an off-shell Q* with virtuality s = Pq (Fig. [I]), pg is the 4-momentum of the 
heavy quark Q, and A4q is the amplitude for producing the heavy quark Q with the same 
3- momentum pq. Here I only present the results for b — > (be) in the S-wave states M: 

Tz(\ — z) 2 r 

D- b ^ So) (z,p ) = N l (1 _ J [6 - 18(1 - 2r)z + (21 - 74r + 68r V 

-2(1 - r)(6 - 19r + 18r 2 )^ 3 + 3(1 - r) 2 (l - 2r + 2r 2 )z A ] (2) 
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for the ground state denoted by B c , and 
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6^bc( 3 Si) 



3N- 



rz(l — z)' 



2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )z 2 



;i - (l-r)z 

-2(1 - r)(4- r + 2r 2 )z 3 + (1 - r) 2 (3 - 2r + 2r 2 )^ 4 



(3) 



for the first excited state denoted by B*, where N = 2a s (2m c ) 2 |-R(0)| 2 /(8l7rm 3 ) and r = 
m c/( m b + m c)- The results for longitudinally and transversely polarized B* mesons can be 
found in Ref. |i~6| , and the results for P-wave can be found in Refs. [[L9|,[2(J. The fragmentation 
functions for charm quark into charmonium can be obtained by putting r = 1/2 in the 
above expressions, and can be found in Ref. |3| JT9|fil| . Gluon fragmentation functions into 
charmonium and bottomium can be calculated in a similar fashion, and the results are 
in Ref. for S-wave, Ref. |22] for P-wave, and Ref. [23|| for l D2 (see also Refs. [24| for a 
derivation from a field theoretical definition). Likewise, the photon and lepton fragmentation 



functions into charmonium were calculated in Ref. 25 



A couple of remarks about these PQCD fragmentation functions are in order, (i) The 
scale of the PQCD fragmentation functions calculated should be of order of the heavy quark 
mass rriQ. We choose it to be mq + 2m q , which is the minimum virtuality of the fragmenting 
quark. In the case of gluon fragmentation functions, the scale is set at 3mg. These fragmen- 
tation functions obey the usual Altarelli-Parisi evolution equations such that fragmentation 
functions at higher scale can be obtained by evolving the Altarelli-Parisi equations, (ii) The 
inputs to these PQCD fragmentation functions are simply the masses of charm and bottom 
quarks, and the wavefunctions at the origin (also derivatives of the wavefunctions and color- 
octet wavefunctions for higher orbital states). All of these quantities can be reliably obtained 
from potential models or from lattice simulations. Since the inputs can be reliably obtained 
from other sources, these PQCD fragmentation functions have high predictive power. 

In Sec. II the fragmentation approach for calculating the B c meson production is de- 
scribed, and the same approach has been used to calculate the fragmentation contribution 
to the ip and if)' production, as will be summarized in Sec. III. Sec. IV describes a fragmen- 
tation model. 



II. B c Meson Production 

The fragmentation approach has been used in calculating the fragmentation contribution to 
the production of be, charmonium, and bottomium. This approach is based on factorization, 
in which the production process is separated into the production of high energy partons 
(quarks and gluons) and the fragmentation of these partons into the meson. In this section, 
I summarize the results from a series of studies on the production rates of the S-wave 

and P-wave (bc) mesons, as well as, the inclusive production of B c . (bc) mesons belong 
to another heavy-quark bound state family, which is made up of a b antiquark and a c 
quark. The spectroscopy for the spin-orbital states is similar to that of charmonium and 
bottomium, and (be) can be obtained by interpolating between charmonium and bottomium. 



According to potential models p6[, the IS, IP, ID, 2S, and possibly the whole set of 2P states 



lie below the BD threshold. A very peculiar feature of the excited (bc) states is that the 
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annihilation channel is suppressed relative to the electromagnetic or hadronic transitions 
into lower-lying states. Therefore, when an excited state is produced it will cascade into 
the ground state with emission of photons and pions. Hence, all the IS, IP, ID, 2S, and 
probably 2P states contribute to the inclusive production of the ground state B c meson. The 
D-wave fragmentation functions are not available yet but they are expected to contribute 
only a very small fraction. 

The differential cross section for producing a (be) meson in a spin-orbital state H is given 

by 



da 
dpr 



(pp^H(p T )X) = Yl J dx 1 dx 2 dzf i/v {x 1:l i)f j/p {x 2 ^) 



da 

xD- b ^ H (z,fi) + -> g(pr/z)X,fj,) D g ^ H (z,n) 

apx 



-7— b(pr/z)X, /i) 
dpr 

(4) 



where /j/ p (x)'s are the parton distribution functions, da's are the subprocess cross sections, 
and D^h{z, /u)'s represent the parton fragmentation functions at the scale fi. The factor- 
ization scale /i is chosen in the order of the pt of the parton, so as to avoid large logarithms 
in da, and the resulting large logarithms of order [i/m b in D^h{z, /i) can be summed by 
the Altarelli-Parisi equations. The gluon fragmentation functions at the initial scale are 
a s suppressed relative to the b fragmentation functions, so we simply take the initial gluon 
fragmentation functions to be zero. This is justified since the majority of the gluon frag- 
mentation comes from the Altarelli-Parisi evolution, and we called it the induced gluon 
fragmentation functions ||. The b and gluon fragmentation functions into a (be) state H 
satisfy the following evolution equations 

d dy /"i dy 

V^-D- b ^ H (z,n) = / —P- b _ ¥ - b (z/y,fi)D- b _ ¥H (y,n)+ —P b ^ g (z/y,fi)D g _ H (y,n), (5) 
o[i j z y J z y 

/"i dy /"i dy 

V^-D g ^ H (z,ii) = —P b (z/y,fi)D b ^ H (y,fi)+ —P g ^ g (z/y, fi) D g ^ H (y, fi) , (6) 
on Jz y Jz y 

where P%^j are the usual Altarelli-Parisi splitting functions. The initial conditions to the 
above equations are simply the initial b fragmentation functions and the initial gluon frag- 
mentation functions, which are set to zero. We can also examine the relative importance of 
these fragmentation functions. The initial D b ^ H (z, /io) is of order a 2 s , even when it is evolved 
to a higher scale it is still of order ct 2 s . In contrast, the initial D g ^ H (z, /i ) is of order a^, but 
at a higher scale /1 it is of order ln(/i/yu ) via the Altarelli-Parisi evolution. Therefore, at a 
sufficiently large scale the induced D 9 ^h(z, /i) is as important as the b quark fragmentation. 

The resulting px spectra for the S-wave and P-wave states are shown in Fig. |2|. We can 
now predict the inclusive production rate of the ground state B c meson. We add up the cross 
sections from all S-wave and P-wave production, and the inclusive cross sections of B c as a 
function of p^ m (B c ) is shown in Table | at the Tevatron. Variations with factorization scale 



between fi R / '2 to 2fi R , where /ir = \jp\ (parton) + ml, are also illustrated. The variation is 
at worst a factor of two, and substantially reduced at p™ m > 10 GeV. With a production 
cross section of about 5 nb and 100 pb _1 integrated luminosity at the Tevatron there are 
about 5 x 10 5 mesons. The detection mode for the B c meson will be B c ^ J/ip + X, 
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where X can be a tt + , p + , or £ + vi, and J/ip can be detected easily through its leptonic 
decay modes. When X is e + z/ e or p^v^, it will be a striking signature of three-charged 
leptons coming from a common secondary vertex. The combined branching ratio of B c — > 
J/ip£ + vi — > £' + £'~£ + V£ (£,£' = e, p) is about 0.2%. This implies that there will be of order 
10 3 such distinct events for 100 pb _1 luminosity at the Tevatron. However, this mode does 
not afford the full reconstruction of the B c . If X is some hadronic states, e.g., pions, the 
events can be fully reconstructed and the B c meson mass can be measured. The process 
B c — ► J/-?/; + 7r + — ► £ + £~n + is likely to be the discovery mode for B c . Its combined branching 
ratio is about 0.03%, which implies about 300 such distinct events at the Tevatron with a 
luminosity of 100 pb _1 . The production rates of B c at the LHC will be of order 10 9 with 
100 fb -1 luminosities, promising a very exciting experimental program on B c mesons at the 
LHC. 

There also exist complete 0(a^) calculations on the production of S-wave B c mesons 
27|,|28|. There is a controversy whether the fragmentation contribution will dominate the 



Table I: The inclusive production cross sections for the B c meson at the Tevatron including the 
contributions from all the S-wave and P-wave states below the BD threshold as a function of 
prp m (B c ). The acceptance cuts are pt(B c ) > 6 GeV and |y(-B c )| < 1- 
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0.87 


1.16 
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0.26 
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0.26 
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0.098 
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production of B c mesons at the large px regions over the non-fragmentation (recombina- 
tion) contribution. The controversy arises because the set of fragmentation diagrams is a 
gauge-invariant subset of the whole set of Feynman diagrams at the order a A s . So there is 
a competition between the fragmentation diagrams and the recombination diagrams. Nev- 
ertheless, the bottom line is that the fragmentation approach identifies the correct scale for 
the fragmentation diagrams, and should give a lower bound on the production cross section 
of B c mesons. 



III. J/ip and ip' Production 

In this section, I highlight some recent activities on J / ip and ip' production at the Tevatron. 
Before the recent CDF data 0, the dominant mechanism of J/if> production at high px 
region was believed to be the fusion mechanism, gg — > Xcj9, followed by the radiative decay 
°f Xcj J /ip + 1) while for ip' production the dominant mechanism is gg — > ip'g because 
of the absence of Xcj(2-P) states. But the CDF measurements on J/ip and ip' exposed 
large discrepancies between theoretical predictions and experimental data. The discrepancies 
demonstrate that there must be either some other unknown production mechanisms or simply 
that perturbative QCD is not valid in this case. 

In order that perturbative QCD is still the means to understand the production of heavy 
quarkonia, it is advantageous to consider higher order contributions, which are more impor- 
tant than the lowest order fusion mechanism. It was shown explicitly in Ref. || that the 
contributions to J/ip and if)' production from gluon, charm quark, and photon fragmenta- 
tion are more important than the lowest order fusion mechanisms beyond certain values of 
Pt (see Fig. |3|.) Among all the fragmentation contributions that are relevant to J/ip pro- 
duction, the largest one comes from gluon fragmentation into Xcj followed by the radiative 



decay Xcj — > J/V + 7 [p2| . The gluon fragmentation D g _> XcJ (z) consists of two pieces, one 
of which is the color-singlet part of order o? s and the other piece is the color-octet part of 
order a s . When the fragmentation contributions are included, the theoretical prediction 
matches the experimental data within a factor of 2 - 3 |8|,p|,p!0| (see Fig. which is within 
the uncertainties from the mass of charm quark, the factorization scale, higher order QCD 
corrections, and relativistic corrections. 

While anomalous J/ip production seems to be solved, however, the data for ip' produc- 
tion is still a factor of 20 - 30 above the theoretical prediction, even after including the 
fragmentation contributions (see Fig. [|). The %cj(2P) states are predicted to be above the 
DD threshold and therefore do not contribute to ip' production. This discrepancy is some- 
times referred as the ip' anomaly. Of course, there have been speculative solutions to the 
anomaly. The most obvious solution is the hypothesis that Xcj(2-P) states are metastable 



such that they decay with appreciable branching ratios into ip' PH|JI2]JI3]JI^| . According 



to potential models, the Xcj(^P) states are above the DD threshold, but the decay of 
Xcj —* DD might be suppressed due to a D-wave suppression. Therefore, an appreciable 
fraction of %cj(2P) can decay into ip'. In order to explain the ip' anomaly, a branching ratio 
B(Xcj(2P) ~ > ip' + 7) ~ 5 — 10% is needed. However, such a large branching ratio is unfla- 
vored by potential models. There is also another mechanism due to Braaten and Fleming 
who proposed that the production is via gluon fragmentation into a color-octet 3 S[ 
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Figure 3: Contributions to the differential cross section for inclusive J /if) production at the 
Tevatron: fragmentation into if) (solid curves), and the leading order contributions (dashed curves). 
(Taken out from Ref. ||.) 



state, which then nonperturbatively emits a pair of gluons to make the transition into if)' 
state. This mechanism is suppressed by powers of v, which is the relative velocity between c 
and c inside the charmonium. But on the other hand, this fragmentation is flavored by two 
powers of a s compared to the color-singlet fragmentation function (af). It means that this 
fragmentation mechansim could be potentially large because the corresponding fragmenta- 
tion function is only of order a s , though suppressed by powers of v. The major uncertainty 
is the determination of the nonperturbative parameter associated with the soft emission of 
the gluon-pair in color-octet + state. This mechanism can be tested rather easily because 



the if) 1 produced will be entirely transversely polarized [115 , 121 , and the polarization can be 



easily measured experimentally by looking at the angular distribution of the muon pair in the 
decay of if)'. If this mechanism is tested to be important, it will significantly affect all hadro-, 
photo-, and electro- production of charmonium. In fact, a very recent analysis indicated that 
only about one third of the prompt J / if) comes from Xcj decays, while the rest is from direct 
J I ip production. It means that there is another important production mechanism other than 
the gluon fragmentation into Xcj's followed by the radiative decay of Xcj's into J /if). This 



might be the hint showing the importance of Braaten-Fleming's color-octet mechanism []T5 
in J / if) production as well. 



IV. A Fragmentation Model 



In this section I describe a fragmentation model [16,^,18] for the fragmentation of heavy 
quarks into heavy-light mesons. This model is based on the PQCD fragmentation functions 
that are presented in Eqs. (|2|) and (|3|) for S-wave mesons, and in Ref. |19| for P-wave mesons. 
But since most of the experimental data are on S-wave states, I concentrate on Eqs. (0)-(@). 
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Figure 4: Preliminary CDF data for prompt ip and ij)' production with theoretical predictions of 
the total fragmentation contributions (solid) and the leading order fusion contributions (dash). 



Let me rewrite the fragmentation functions here: 



D 



Q^QqQ-So) 



>,/*o) 



N- 



rz(l — z) 



l-r)z) 



6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )z 2 



2(1 - r)(6 - 19r + 18r 2 )z 3 + 3(1 - r)\l - 2r + 2r 2 )z 4 



(7) 



Dq^ss^z, lio) = 3N _ V _ > [2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )z 2 

-2(1 - r)(4 - r + 2r 2 )z 3 + (1 - r) 2 (3 - 2r + 2r 2 )z 4 

where iV, instead of being given in terms of the wavefunction and the running coupling 
constant, is now treated as a free parameter governing the overall normalization together 
with another free parameter r, r being the mass ratio m q /(mQ + m q ). This model, in the 
limit r — > 0, approaches the fragmentation of a very heavy quark into a heavy-light meson. 
In reality, the heavy quark is either b or c, and the light component is u, d, or s, and 
therefore r is a small parameter. In principle, when the lighter constituent quark becomes 
a light quark, there is that nonperturbative physics involved in the fragmentation process. 
But we do expect our PQCD fragmentation functions with iV and r as free parameters can 
at least provide a qualitative picture and hence a reasonable model for fragmentation into 
heavy-light mesons. This model is suitable for b — > B,B*,B**, ... and c — > D,D*,D**, ... 
mesons. 

This fragmentation model has certain advantages over previous models |29| in the liter- 



ature. First of all, this model lies on a firm basis of PQCD. It is rigorously correct in the 
limit when m q 3> Aqcd and higher order corrections can be systematically calculated. The 
spirit of our model is the continuation of m q to a value close to Aq CD . The most obvious 
advantage of this model is the ability to predict different results for different spin-orbital 
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states with only two parameters, in contrast to the Peterson model. Another advantage is 
that the fragmentation functions for the same orbital angular momentum share the same 
paramter N, as shown in Eqs (^) - (Q). This is a substantial improvement when ratios of the 
fragmentation functions are measured, in which the N dependence cancels out. For example, 
our model can predict the ratio Py = V/(V + P) as a function of r only, where V is the 
vector meson and P is the pseudoscalar. 

Another theoretical issue is that the PQCD fragmentation functions are consistent with 
heavy quark symmetry, which I explain in more detail next. According to an analysis using 
the heavy quark effective theory (HQET) []3D| , the fragmentation function for a heavy quark 
Q into a hadron Hq containing a single heavy quark Q at the heavy quark mass scale mg is 

D Q ^ HQ (z) = ^a(y) + b(y) + 0(r) (9) 

which is a heavy quark mass expansion in r = — — — - and y is a rescaled variable of z, 

y = 1 ~ ( ',*~ r ' >2 • The leading term is of order 1/r, i.e., mg, while the next-to-leading term is of 
order r°. The PQCD fragmentation functions in Eqs. (^) and @ can be expanded in powers 
of r and reexpressed in terms of y, as 



D Q ^ S0) (Z) = M^(l (3y 2 + % + 8) + (3y 3 + 15y2 + 8y _ 8) + .. 

D Q ^ Qq{3Sl) (z) = ^H!(f(3 2/ 2 + % + 8)-3( 2 / 3 + 2 / 2 -8y + 8) + ...) 



(10) 



The above expansions in powers of r are consistent with the HQET analysis |30| that the 
leading term is order 1/r. Actually, if examined more carefully, the leading terms of the 1 So 
and 3 Si expressions are exactly in the ratio of 1:3, which is the value predicted by heavy quark 
spin symmetry. The next-to- leading terms are not in the ratio of 1:3, and they explicitly 
break the spin symmetry. This fact also prompted us to derive the PQCD fragmentation 
functions independently from the HQET Lagrangian fll8| . By using the leading and the 
next-to-leading terms of the HQET Lagrangian, we succeeded in obtaining the same results 
as the heavy quark mass expansions in Eqs. (0). Therefore, with the consistency with heavy 
quark symmetry and HQET we have more confidence in applying our PQCD fragmentation 
functions as a fragmentation model for heavy quark fragmenting into heavy-light mesons, 
namely, b — > B,B*,B** and c — > D,D*,D**,... When more data on P-wave mesons are 
available, comparisons with the P-wave fragmentation functions can also be made. For the 
moment the data are more or less entirely on the S-wave states. I shall demonstrate a couple 
of comparisons between the predictions by this model and the experimental data. We shall 
look at Py and (z). 

(A) Py for the charm system is defined as Py = D*/(D + D*), which is a measure of the 
population of D* in a sample of D and D* mesons. Since fragmentation is the dominant 
production mechanism for charm mesons and D, D* mesons dominate, Py can be expressed 
in terms of the fragmentation functions 

J dzD c ^D*(z) 

V= fdzD c ^ D (z)+fdzD c ^ D .(z) ' (U) 

which is a function of r only. The prediction by the PQCD fragmentation model is shown 
by the solid curve in Fig. |5|. At r = 0, the heavy quark mass limit, Py = 0.75, which is 
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> 




Figure 5: The ratio Py = V/(V + P), where V denotes the vector state (D*,B*) and P denotes 
the pseudoscalar (D,B). The data for D,D* system and for B,B* system are shown. 



exactly the value given by the naive spin counting. At r > Py is always smaller than 
0.75, which implies that D* is produced less than given by heavy quark spin symmetry. This 
can be understood in terms of the mass splitting between D and D* mesons, which can be 
accounted for by the a^G^/M term in the HQET. 



A compilation of data on Py can be found in Ref. [3T|, in which the updated branching 
ratio B(D* + -> D°tt + ) = 0.681 ± 0.016 was used, and the average Py = 0.646 ± 0.049. This 
value of Py also indicates that D* mesons are produced less than it should be as given by 
heavy quark spin symmetry. For the charm quark we take m c = 1.5 GeV and the light 
constituent quark mass to be 0.3 GeV, therefore r = m ' lght — = 0.167. The data is then 
plotted on the graph and very good agreement is obtained. Recently, the data for B, B* 
system has also been available with P V (B, B*) = 0.76 ± 0.08 ± 0.06 @. We take m b = 4.9 



GeV and m light = 0.3 GeV again, therefore r = 0.06. The prediction is still less than la 
from the data point. From the figure we can see that if we choose a smaller value of mnght, 
say 0.2 GeV, we could even get a better agreement. The errors in Py certainly allow us to 
vary mug^t more than 0.15 GeV such that the prediction is still within la. Or, we can use 
the experimental value of Py to fix the parameter r. 

B) (z) is the average longitudinal momentum fraction that is transferred from the heavy 
quark to the meson. In terms of fragmentation functions, (z)^ at a scale \x is given by 

,\n _ JdzzD c ^ D *(z,fj) 
J dz D c ^ D *(z,ii) 

Experimentally, the inclusive c — > D* channel was measured at LEP, at CLEO, and at 
ARGUS. The D * given in Eqn. (|i~2"D is the ratio of the second to the first moments of 
the fragmentation function at the scale \x. Since the anomalous dimensions of the moments 
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Figure 6: The average (z)^ for c — > D* and for b — > .B* fragmentation versus the scale The 
experimental measurements from LEP (/i = m^/2) and from CLEO/ARGUS (// = 5.3 GeV) are 
shown. 

are known explicitly, the scaling behavior of can be determined to be 

where 7 = — 4CV/3, CV = 4/3, 6 = (lliV c — 2n/)/3, iV c = 3, n/ is the number of active 
flavors at the scale /x, and (z)^ is the value determined at the initial scale fiQ. Taking the 
inputs: m c = 1.5 GeV, m Uid = 0.3 GeV, /i = fn c + 2m Utd = 2.1 GeV, we have r = 0.167, and 

(z)c^d* = and (z)^Z^* = 0.50. The variation of (z)c_> D * and {z)%->b* as functions of 
/i are shown in Fig. [?], where we chose m b = 4.9 GeV. The curve for (z)^ B * was also shown 
because we want to demonstrate that using the same m Ujd and m b = 4.9 GeV, the results 
predicted also agree with the data for the bottom quark fragmentation at LEP. 

The measured quantity is (xe), which is the energy of the meson relative to one half 
of the center-of-mass energy of the machine, and xe should be a good approximation to z. 
At LEP, the average value of (xe)c->d* = 0.504 ± 0.0133 p3] . For the bottom quark, only 
the inclusive hadron production has been measured. But we expect that {xe)i)->b* should 
be close to (xE)b^H b , where H b is a bottom hadron, because the b — » B* is the dominant 
fragmentation mode of the bottom quark. The average value of (xE)b^H b = 0.694 ± 0.0166 
34|. Also, we have data on c — > D* from CLEO and ARGUS Combining the CLEO and 



ARGUS data we have (xe) c ^d* = 0.648 ± 0.043. The scale of the measurements is taken to 
be one half of the center-of-mass energy of the machines, so it is m^/2 at LEP and 5.3 GeV 
at CLEO/ARGUS. These data are shown in Fig. || Excellent agreement is demonstrated. 
The only inputs to these comparisons are simply r and hq. Once they are fixed, (z) Mo can be 
calculated by Eqn. ([12]) and evolved by Eqn. (|i~3"|) to any scale fi. The results show agreement 
at two different scales. This is a big contrast to Peterson fragmentation model, which fits to 
different values of mass ratio eg at different scales. 
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V. Conclusions 



In this proceedings I have summarized some recent work on PQCD fragmentation functions. 
The ip production at the Tevatron is the first evidence showing the importance of parton 
fragmentation in quarkonium production in the large px region. I have also shown the re- 
sults of the production of be mesons by the fragmentation approach. Finally, a fragmentation 
model based on the PQCD fragmentation functions is advocated to describe the fragmenta- 
tion of heavy quarks into heavy-light mesons. This model lies on a firm basis of PQCD and is 
consistent with HQET, and was successfully applied to fit data on Py and (z)^. Other work 
includes the estimation of the strange-quark mass parameter m s in the B s and B* system, 
in which the probability of a b quark going into a stranged B meson is fitted to the value of 
m s [36], and a value of about 300 MeV was obtained. There was also a calculation [37] of 



the Falk-Peskin variable w 3 / 2 [0 using the PQCD fragmentation functions. Another work 
was the extension to the fragmentation of a heavy quark into a baryon containing two heavy 
quarks 
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